Non-Markovian Dynamics of Quantum Discord 
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We evaluate the quantum discord dynamics of two qubits in independent and common non- 
Markovian environments. We compare the dynamics of entanglement with that of quantum discord. 
For independent reservoirs the quantum discord vanishes only at discrete instants whereas the 
entanglement can disappear during a finite time interval. For a common reservoir quantum discord 
and entanglement can behave very differently with sudden birth of the former but not of the latter. 
Furthermore, in this case the quantum discord dynamics presents sudden changes in the derivative 
of its time evolution which is evidenced by the presence of kinks in its behavior at discrete instants 
of time. 

PACS numbers: 03.65.Ud, 03.65.Yz, 03.67.Mn 



I. INTRODUCTION 

Entanglement is a kind of quantum correlation (QC) 
that has been playing a central role in quantum informa- 
tion and communication theory pQ. However there are 
other nonclassical correlations apart from entanglement 
[2HI] that can be of great importance to these fields. In 
order to characterize all nonclassical correlations, Ollivier 
and Zurek introduced what they called quantum discord 
[2] . This measure of quantum correlations captures a fun- 
damental feature of classical bipartite states - when the 
discord is zero the information is locally accessible and 
can be obtained by distant independent observers with- 
out perturbing the bipartite state. Although a vast lit- 
erature exists on the study of entanglement just recently 
the other quantum correlations received due attention 
[T2"] . A motivation for the study of these correlations, for 
example, is the recent discovery that nonclassical correla- 
tions other than entanglement can be responsible for the 
quantum computational efficiency of deterministic quan- 
tum computation with one pure qubit (DQC1) In 
this context, the quantum discord could be a new re- 
source for quantum computation. 

However, realistic quantum systems are not closed and 
therefore it is of fundamental importance to study the 
quantum correlations when the system loses its coher- 
ence due to interactions with the environment |13j . The 
entanglement dynamics in open quantum systems was 
broadly studied in the literature but not much exists 
about the effect of the environment on quantum discord 
[T01 - H2] . A peculiar aspect of the entanglement dynamics 
is the well-known "entanglement sudden death" (ESD) 
phenomenon [2J [T3] . This process described the finite- 
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time disentanglement of two parts that interact with ei- 
ther independent [15H17) or common environments |18l - 
[20] . In a previous work [TU] we observed that, even at 
finite temperatures, under a dissipative Markovian pro- 
cess, the quantum discord is immune to "sudden death" . 

Despite the term "entanglement sudden death" sounds 
mysterious, it is important to note that there is no 
dynamical distinction between separable and entangled 
states, since the quantum states can, in general, evolve 
back and forth across the boundary between distinct full- 
dimensional subsets of the space of the density matrices 
which contain separable and entangled states. 

The studies on the entanglement dynamics, that ini- 
tially were restricted to Markovian approximations, have 
recently been extended to consider non-Markovian en- 
vironments [TTl [T9H2TJ . In this case, given the mem- 
ory stored in the environment, some of the initial en- 
tanglement that is lost during the dissipative dynamics 
can return to the qubits. This phenomenon is known 
as "sudden-birth of entanglement" (SBE), which in the 
light of what has been said in the previous paragraph 
should not present much of a surprise to us. Nevertheless 
a question still remains; what happens to the quantum 
discord in this situation? Since quantum discord exists 
even without entanglement, does it present sudden death 
or even sudden birth? In order to answer these questions 
we study the quantum discord of two qubits coupled to 
non-Markovian dissipative environments. 

In this paper we evaluate the quantum discord dynam- 
ics for a dissipative non-Markovian process. For inde- 
pendent environments, when the qubits are subject to 
amplitude damping, we show that it only vanishes at 
discrete instants of time, each within the time interval 
when the reduced quantum state becomes pure and sepa- 
rable, and, consequently, the entanglement vanishes. For 
a common reservoir the quantum discord behavior can 
be very different from that of the entanglement. While 
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the entanglement dynamics presents damped oscillations 
with or without sudden death, the quantum discord is 
almost always positive and presents isolated kinks (cus- 
pids) at which there is a jump in its derivatives. The 
latter behavior is at clear variance with what happens to 
the entanglement. 



II. QUANTUM DISCORD 

Entanglement is not the only measure of quantum cor- 
relations and therefore an interesting approach was intro- 
duced in [31 [3J to attempt to quantify all the nonclassical 
correlations present in a system besides entanglement. 
The defined quantity - the quantum discord - is given by 
the difference between two expressions of mutual infor- 
mation (MI) extended from classical to quantum system. 

The total correlation between two classical systems A 
and B, whose state is described by a joint probability dis- 
tribution p (A, B) , can be obtained by a measure of the 
MI, X(A : B) = H (A) + H (B) - H (A, B), where H (•) 
denotes the Shannon entropy H (p) = — J^jkPjk ^°&2 Pjk 
PP. This classical MI can be rewritten as the equivalent 
expression J (A : B) = H (A) - H (A\B) through the 
Bayes rule [24] , where the conditional entropy H (A \ B) 
quantifies the ignorance about the state of A when one 
knows the state of B. For a quantum system repre- 
sented by a bipartite density operator p, the Shannon 
entropy functional is replaced by the von Neumann en- 
tropy, S (p) = — Tr (plog 2 p), which is the first quantum 
extension of the classical ML We denote it I (p) . 

Another route to generalizing the classical MI to the 
quantum case is to use a measurement-based condi- 
tional density operator [2|. If we restrict ourselves 
to projective measurements performed locally only on 
system B described by a complete set of orthogonal 
projectors, {11/,.}, corresponding to outcomes fc, the 
quantum state after a measurement changes to pk — 
[(I ® n k ) p (I ® II fe )] /Tr (I ® H fe ) p (I ® n fc ) , where I is 
the identity operator for system A. With this condi- 
tional density operator, a quantum analogue of the con- 
ditional entropy can then be defined as S (p | {ITfc}) = 
SfcPfe^ (Pfc)> an d the second quantum extension of the 
classical MI may be found, J (p | {LT fc }) = S (p- 4 ) — 
S (p | {ilfc}). Projective measurements on system B re- 
move all nonclassical correlations between A and B, but 
the value of J (p | {LI fe }) depends on the choice of {II fe }. 
Therefore, to ensure that it captures all classical correla- 
tions, we need to maximize J over all {Lt^}. This quan- 
tity, Q(p) = su P{n fe } J {P I {Hfc})j is interpreted explic- 
itly by Henderson and Vedral [5] , as a measure of classical 
correlations. The quantum discord is then defined as 



D{p)=X(p)-Q{p), 



(1) 



and provide us with information on the quantum nature 
of the correlations between two systems, such that it is 
zero only for states with classical correlations [21 [3J and 



nonzero for states with quantum correlations. Although 
quantum discord is equal to the entanglement of forma- 
tion for pure states , it is not true for mixed states, since 
some states present finite quantum discord even without 
entanglement [3J. 

It is important to note that to calculate the classical 
correlations one can consider arbitrary POVM measure- 
ments as Henderson and Vedral did in [3J. However, for 
two qubits, which is our case, Hamieh et al. [55] show 
that the projective measurement is the POVM which 
maximizes the classical correlations. 



A. Analytical expression for quantum discord 

To evaluate the quantum discord dynamics presented 
in this article we determine an analytical expression for 
a subclass of the X structured density operator. We con- 
sider a density matrix as given by 



P(t) 



a w 

b z 

z b 

w d 



(2) 



where the coherences are real numbers and the element 
P22 = p3S- It is easy to see that for this expression of 
p(t) the condition S{p A ) = S(p B ) is satisfied and there- 
fore the measurement of classical correlations assumes 
equal values, irrespective of whether the measurement is 
performed on the subsystem A or B 3J. To reduce the 
difficulty to compute the quantum discord we need to be 
able to maximize the classical correlation Q{p). This can 
be done analytically if one notes that a general one-qubit 
projector can be written as a function of two angles, since 



Q(p) = S(p A )-F(0,< 



(3) 



where 



F(fi, 



inf 



fe=l,2 



k \ u i 1 



k ' 



Pk{6A) 



(4) 

with p k {e,cj>) = Tr{nf(M)^ B nf(M)} and the pro- 
jectors nf(0,0) = I® \k) (k\ (k = 1,2) defined by the 
orthogonal states: 



|1) = cos 6 |t> +e** sin 9 ||) 
|2) = sin 6 |t) -e i0 cos 9 \i) 



(5) 



We begin noting some peculiar properties of F(9, <j>) 
when p AB is given by Eq. (2]). Given the structure of the 
density matrix the critical points of F(9, </>), i.e. the set of 



values of 9 and 6 such that 



dF(6>,0) 



and 



dF(9,, 



0. 



80 ~ u d<f> 

do not depend on the elements of the density matrix. 
For 9 = tit, with n € Z we have a set of critical points 
and in this case the function F{9, 4>) does not depend on 
the angle ip. Another set is given by 9 = rrij and <f> — 
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with m, n € Z. Thus, with this observation, using 
the quantum version of the mutual information I and 
Eq. ([I]), it is straightforward to compute an analytical 
expression for the quantum discord: 



D{p) = min{L>i,L> 2 } 



(6) 



where 

D t = S(p*)-Sir A *)-«^2 



■4, C,...^, ....... / " \ _ 61og ' 



and 



a + b 



a + b 



(7) 



^2 = SVVS(p^)-A+log 2 A + -A_log 2 A_,(8) 

with A± = \ (1 ± T) and r 2 = (a - df + A(\z\ + \w\) 2 , 
which has been numerically verified for any density op- 
erator with the same structure as in Eq. ([2]). 



III. EXACT DISSIPATIVE DYNAMICS OF 
QUANTUM DISCORD 

In this article we study a system whose dynamics is 
described by the well-known damped Jaynes-Cummings 
model. We consider two distinct situations, the inde- 
pendent and common environment. In the former, each 
qubit is coupled to its own reservoir since the dissipative 
processes occur independently. In the common environ- 
ment case, on the other hand, we consider only one bath 
for both qubits. We suppose that the qubits are cou- 
pled to a single cavity mode which in turn is coupled 
to a non-Markovian environment that initially is in the 
vacuum state. In this case we can say that the inter- 
action Hamiltonian reduces the amplitude of motion of 
each qubit state, and this justifies why such a process is 
known in the literature as the amplitude damping chan- 
nel. The solutions of these simple models have recently 
been used to study the non-Markovian effects on the dy- 
namics of entanglement jTTJ, [5D] . The environments are 
represented by a bath of harmonic oscillators, and the 
spectral density is of the form 



J(w) 



f 



7oA 2 



2tt (w ~ + A 2 ' 



(9) 



where A is connected to the reservoir correlation time 
tb by the relation tb w I/A, and 70 is related to the 
time scale tr over which the state of the system changes, 
tji w l/7o- Here we will consider the strong coupling 
limit, i. e. tr < 2tb- 

For independent amplitude damping channels the two- 
qubit Hamiltonian can be written as 



tt _ , V, ,(*)„(*)*„(») 

H - uj a + cr_ +2^ w fe a k a k 



where i?W — Y^kQk a k with 9k being the coupling 
constant, Wq is the transition frequency of the i-th qubit, 
and <T^j!' are the system raising and lowering operators 
of the i-th qubit. Here the index k labels the reservoir 

field modes with frequencies oJf 1 , and a£ (c4 ) is their 
usual creation (annihilation) operator. Here and in the 
following the Einstein convention sum is adopted. 

For common environments, on the other hand, we have 
that the system raising (lowering) operator of each qubit 
is coupled to the same environment operator B (B*). In 
this case we have one bath coupled to both qubits, and 
the Hamiltonian is given by 

H = J£ ) afa^+£u k a k J !a k +(afB + ( T^B^ . (11) 

k 

We have considered two identical atoms equally cou- 
pled to the reservoir. In this case, the dynamics of the 
two qubits occur in two completely decoupled subspaces, 
generated by {|00> , |+) = (|10) + j01>) /-y/2, |11)} and 
(|-) = (| 10) - |01>) /V2}. Using this fact, Mazzola et 
al. [20] connect the problem with a three-level ladder 
system [35] and through the pseudomode approach |26j . 
the dynamics can be evaluated without any approxima- 
tions. 

The initial state considered in this paper is the Bell-like 
state: 



1 00) + \j\ - a 2 |11) 



(12) 



and thus the density matrix of the atomic system has the 
form of Eq. ([2| with the dynamics of its matrix elements 
given in [T7j for independent reservoir and (2Ql [23] for 
common reservoir. 

In Fig. we plot the quantum discord as a function 
of the scaled time 70^ in the strong coupling regime, A = 
O.I70 for the common reservoir case and A = O.OI70 for 
independent reservoirs. The initial state used in Fig. ([T]) 
is given by Eq. (12) with a 2 — 1/3. The analytical 
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solution of the quantum discord is the minimum value 
assumed by the functions D\ (squares) and D 2 (circles) , 
given by Eq. |7| and Eq. Q respectively. In Fig. it 
is represented by the black (solid) line. 

We have observed that, for the values of •jot indicated 
by the doted lines in Fig ([!]), the angles that minimize 
the discord change, point out to a "sudden change" of 
discord as previously noted in |llj for Markovian envi- 
ronments. In the region between two consecutive dotted 
lines the angles remain unchanged. As shown in Fig. lb, 
the same analysis applies to independent reservoirs, but 
no "sudden change" is observed in this case. Here it 
should be emphasized that what we are calling "sudden 
change of discord" is actually a signature of a jump of 
the time derivatives of that function at specific instants. 

In order to compare the discord dynamics with the 
entanglement dynamics we used the entanglement of for- 
mation (EoF) [57] as a measure of entanglement. For two 
qubits, the EoF dynamics can be written as a function 
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FIG. 1: (Color online) Analytical (Di-squares, and D2- 
circles) and numerical (solid line) dynamics of discord for (a) 
a common reservoir with A = O.I70 and (b) independent reser- 
voirs with A = O.OI70. The dotted lines in (a) indicate the 
value of 70 1 where the sudden change occurs. The two-qubits 
initial state used here is given in Eq. (12 1 with a 2 = 1/3. 



of the concurrence |28j and it is given by 

E(t) = -r(t) io g2 r(* ) - [1 - r(t)} io g2 [1 - r(t)] (13) 

where T(t) = |(1 + y/l - C(i) 2 ) with C(t) being the 
time dependent concurrence. For a density matrix with 
a structure defined as in Eq. ^ we have that C(t) = 
2max{0, Ai(t), A 2 (t)} with Ai(i) = z(t) - y/a(t)d{t) and 
A 2 (i) =w(t)-b(t). 

We begin analyzing the entanglement and quantum 
discord dynamics for independent reservoirs. The en- 
tanglement dynamics, for example, shows different be- 
haviors depending on the initial state of the two-qubit 
system [17]. In the case where a 2 > 1/2 the EoF 
periodically vanishes for the discrete times defined by 
t n = 2 [nn — arctand/A] /d with d = y // 2"/ X — A 2 and n 
an integer [T7]. This behavior is illustrated in Fig. (2c) 
(dashed line) where we observe that the amplitude of os- 
cillation undergoes a decay after each revival. However, 
when a 2 > 1/2 the behavior of the EoF presents two 
different features: (£) there is ESD because the EoF per- 
manently vanishes within finite time intervals [see Fig. 2a 
(dashed line)] and (ii) the revival of entanglement after 
these intervals when the two qubits are fully disentangled 
[see Fig. 2b (dashed line)]. 
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FIG. 2: (Color online) Dynamics of Discord (solid line) and 
EoF (dashed line) as a function of the scaled time for indepen- 
dent reservoirs with A = O.OI 70 a nd the two qubits initially 
prepared in the state of eq. (121 with (a) a 2 — 1/10, (b) 
a 2 = 1/3, and (c) a 2 = 1/2. 



The discord dynamics for all initial entangled states 
(0 < a 2 < 1) is similar to the EoF for a 2 > 1/2, i.e., the 
discord vanishes only at t n when the two-qubits state 
becomes the separable pure state 1 00) [see Fig.2a-c (solid 
line)]. The nonclassical correlations are mediated by the 
reservoir, since there is no interaction between the qubits. 
Furthermore, whereas the entanglement may reappear af- 
ter a time interval within which the EoF is zero, the dis- 
cord is almost always non zero. This result indicates that 
the discord under non-Mar kovian dissipative dynamics, 
likewise the Markovian case , vanishes only at discrete 
instants. This point agrees with the results presented in 
|12) . where the authors show that the states with zero 
QD form a set of measure zero. 

In the case where the two-qubits interact with the 
same environment the entanglement dynamics presents 
two regimes [20]: for a 2 > 1/4 damped oscillations of 
entanglement are observed [see Fig. 3d, dashed line], and 
for a 2 < 1/4 finite time intervals of complete disentangle- 
ment are followed by entanglement revivals [see Fig.3b-c, 
dashed line]. However, while for independent reservoirs 
the discord and entanglement behaviors are similar to 
each other, for common environment they behave very 
differently, as show in Fig.3a-d. For a 2 — 1/2 both the 
discord and entanglement dynamics present the same be- 
havior, but as a 2 decreases the discord exhibits very com- 
plicated damped oscillations with sudden changes more 
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FIG. 3: (Color online) Dynamics of Discord (solid line) and 
EoF (dashed line) as a function of the scaled time for common 
reservoirs with A = 0.170 and two qubits initially prepared in 
the state of eq. {T2I with (a) a 2 = 0, (b) a 2 = 1/10, (c) 
a 2 = 1/5, and (d)c? = 1/2. 
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FIG. 4: (Color online) Dynamics of Discord (solid line) and 
EoF (dashed line) as a function of the scaled time for a 
common reservoir with (a) A = O.I70, (b) A = 70, and (c) 
A = IO70. The two qubits initial state is given by eq. (12 1 
with a 2 = 1/3. 



evident. The difference between these two measures is 
even more drastic when a 2 = 0, where the two-qubit 
state is initially the separable state |11). In this case the 
interaction between the qubits mediated by a common 
reservoir does not lead to the generation of entanglement 
between them [20]. On the other hand, the reservoir- 
mediated interaction leads to the generation of nonclas- 
sical correlations as exposed in Fig. (3a) , showing a "sud- 
den birth" of discord but not of entanglement. Moreover, 
since the reservoir is initially in the vacuum state, when 
the initial state is |00), obtained from Eq. (12) with 



a — 1, no correlation is created since the composite sys- 
tem is in its ground state. It is also interesting to note 
in Fig.3b-c that the decrease of entanglement is accom- 
panied by the increase of discord in some regions. 

In Fig. 4 we plotted the discord (solid line) and the EoF 
(dashed line) as a function of the scaled time j^t for a 
common reservoir with (a) A = O.I70, (b) A = 70, and 
(c) A = IO70. The two-qubits are initially in the state 
(12) with a 2 = 1/3. These results show that the entan- 
glement decay almost exponentially for A = IO70, as in 
the case of Markovian reservoirs. It is expected because 
in this regime A > 270 and therefore the qubit-reservoir 
coupling is weak. In this same weak coupling regime the 
discord has a similar behavior and we note that the num- 
ber of points where the sudden change occur also tends to 
decrease. Besides, as the effective coupling between the 
two qubits is due to the action of the common reservoir, 
the quantum correlations created tends to decrease. 

It is worth mentioning that results very similar to ours, 
for the case of independent environments, have been nu- 
merically studied in [2H] . 



IV. CONCLUSION 

We have studied the quantum discord dynamics of 
two-qubits coupled to common and independent non- 
Markovian environments. We have used the exactly solv- 
able damped Jaynes-Cummings model for zero temper- 
ature environments. We have observed that even when 
the entanglement suddenly disappears and reappears af- 
ter finite time intervals, the quantum discord vanishes 
only at discrete times. For a common environment we 
have observed what is called the sudden change phe- 
nomenon. Actually, the quantum discord between the 
qubits suddenly changes depending on the maximization 
process of the amount of classical correlations between 
them. This fact indicates that this phenomenon could 
be universal which means that, for general initial condi- 
tions and interaction Hamiltonians, we expect that the 
POVM that maximizes the classical correlations would 
abruptly change from one time interval to the other. This 
point will be studied in a future work where more general 
Hamiltonians and initial conditions will be considered. 
Furthermore, we have observed that, in the case of com- 
mon environments, the very different behavior of discord 
and entanglement can arise even for initially separated 
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states. Finally, we have also noticed that, even without 
entanglement, the correlations introduced by the envi- 
ronment are transferred to the two qubits producing a 
finite quantum discord. 
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